Simultaneous cooling of an artificial atom and its neighboring quantum system 
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We propose an approach for cooling both an artificial atom (e.g., a flux qubit) and its neighboring 
quantum system, the latter modeled by either a quantum two-level system or a quantum resonator. 
The flux qubit is cooled by manipulating its states, following an inverse process of state population 
inversion, and then the qubit is switched on to resonantly interact with the neighboring quantum 
system. By repeating these steps, the two subsystems can be simultaneously cooled. Our results 
show that this cooling is robust and effective, irrespective of the chosen quantum systems connected 
to the qubit. 
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Quantum devices using Josephson junctions can be 
used as artificial atoms (AAs) for demonstrating quan- 
tum phenomena at macroscopic scales. Witli states in- 
volving the two lowest energy levels, these devices are 
good candidates for solid-state qubits fH. When using 
their three lowest levels, such a solid-state three-level sys- 
tem, fabricated on a microelectronic chip, can be useful 
for single-photon production Q and lasing Q . 

For single-photon production [§| and AA lasing 
a state population inversion is established for the two 
working energy levels via a third one (i.e., transitions 
|0) |2) |1) in Fig. [1]). Interestingly, the inverse 
process of state population inversion (i.e., transitions 
|1) |2) |0) in Fig.[T]) can be used to increase the oc- 
cupation probability of the ground state and thus lower 
the temperature of the qubit. This idea has been applied 
in a recent experiment [j| to cool a flux qubit. Indeed, 
this is analogous to the optical side-band cooling method 
studied earher (see, e.g., [1| and @). The experiment Q 
shows that the temperature of the flux qubit can be low- 
ered by up to two orders of magnitude with respect to 
its surroundings. This provides an efficient approach for 
preparing a flux qubit in its ground state. 

While the flux qubit was greatly cooled in Q , the noise 
sources surrounding the qubit were not. This is because 
of the weak coupling between the qubit and its environ- 
ment in [3| , where the transition rate between the ground 
and first excited states is small. Below we use a tunable 
AA (to be specific, we choose a flux qubit, but it could be 
another AA) to achieve a strong and switchable coupling 
between the AA and its neighboring quantum system, 
and propose an approach to simultaneously cool both of 
them and not just the AA. Here we consider two typ- 
ical quantum systems to describe the environment sur- 
rounding the AA: (i) a quantum two- level system (TLS), 
which is exactly solvable; and (ii) a quantum resonator. 
Actually, a quantum TLS can describe the noise source 
like a two-level fluctuator and the quantum resonator can 
model the dominant bosons of a thermal bath. In this 



case, the approach is to cool both the fiux qubit and 
such noise sources. This simultaneous cooling of the fiux 
qubit and its neighboring noise sources can significantly 
enhance the quantum coherence of the fiux qubit because 
the cooled qubit is thermally activated very slowly to the 
first excited state, after its neighboring noise sources are 
also cooled. Moreover, the present approach has wide 
applications because the models used here can describe 
other quantum systems. Also, we show that different 
surrounding quantum systems (either a quantum TLS or 
a quantum resonator) give similar results, implying that 
the cooling is robust and effective, irrespective of the cho- 
sen neighboring quantum system. 

Cooling the artificial atom and ground-state 
preparation. — The commonly used fiux qubit 0, S| 
(which is an example for an AA) consists of a su- 
perconducting loop interrupted by three Josephson 
junctions (two equal and one smaller) and pierced by 
a magnetic fiux <I>e. To obtain a tunable AA, the 
smaller junction is here replaced by a SQUID threaded 
by a fiux <I>s [see Fig. [Ila)]. The Hamiltonian can 
be written as 7? = Ppf^Mp -\- Pq/2Mq + U{ifp,ipq), 
with = -ihd/dip, [i = p,q), Mp = 2Cj($o/2^)2, 
and Mq = Mp{l + 47)/4. The potential is U{(pp,ipq) = 
2Ej [1 - cos ifp cos(7r/ -\-^(Pq)]+ 2^Ej [1 - cos(7r/5 ) cos (pq), 
where ipp = {(fi f2)/2, (fiq = {tp-i + <^4)/2, fs = *s/*o, 
and / = ^e/^o + /s/2 (^o is the fiux quantum). To 
drive a resonant transition between states \Ei) and \Ej), 
one can apply a microwave field through the circuit loop: 
$^(i) = $(„"^ cos{uiijt-ir6), with ujij = {Ei-Ej)/h. When 
the microwave field is weak, the time-dependent pertur- 
bation Hamiltonian can be written as V{t) = — /$u,(t), 
where I = —Ic cos ipp sin(7r/ + ^(pq), with Ic = 2nEj/^Q. 
The rate of the state transition between \Ei) and \Ej) is 

Tij (X where tij = {Ei\I^w^\Ej) is the transition 

matrix element. When a neighboring quantum system, 
e.g., a noise source, is coupled to the fiux qubit via a 
fiux variation, then in tij becomes the amplitude of 
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FIG. 1: (Color online) (a) Schematic diagram of an artificial 
atom (A A) produced by a superconducting quantum circuit. 
A symmetric SQUID and two identical Josephson junctions 
with coupling energy Ej and capacitance Cj are placed in a 
superconducting loop pierced by a magnetic flux <l?e (green). 
The two junctions in the SQUID have coupling energy 'yEj 
and capacitance yCj, and the flux (yellow) threading through 
the SQUID loop is Here 7 = 0.5 and Ec = e^/2Cj 
is the single-particle charging energy of the junction, (b) 
and (c): Energy levels of the superconducting A A as a func- 
tion of the reduced magnetic flux / = ^e/^o + fsfi, for 
fs = Os/^o = 0.22 and 0.27, where only the four lowest 
levels are shown and the energy is in units of Ej. (d) and 
(e): Moduli of the transition matrix elements \tij\ (in units of 
Jc^L"') as a function of /, for fs = 0.22 and 0.27. Note that 
each figure in (b)-(e) is symmetric about / = 0.5 and half of 
it is plotted. The vertical (red) dashed lines at / = 0.493 are 
just a guide to the eye. (f) Transition diagram of the AA. At 
nonzero temperatures, the fiux qubit is thermally activated 
from the ground state |0) to the first excited state A 
resonant transition from |1) to the second excited state |2) is 
driven by a microwave field, so as to eliminate the unwanted 
thermal population of and followed by a fast decay to |0). 
While the qubit is cooled to its ground state |0), the A A is 
then switched on, to resonantly interact with a neighboring 
quantum system for a period of time. Repeating these pro- 
cesses, both the qubit and the neighboring quantum system 
can be simultaneously cooled. 



the flux variation. 

In Figs.[Hb)-[IJe), we show the energy levels of the AA 
and the transition matrix elements \tij\ for two values 
of fs- The SQUID gives an effective Josephson coupHng 
energy aEj with a ~ 0.77 and 0.66, respectively. For 
any nonzero temperature, the system will be thermally 
activated from the ground state |0) = \Eo) to the first 
excited state |1) = \Ei). Here we consider the case in 
Fig. [Ijb), with the system working at, e.g., / = 0.493. 
As shown in Fig. 1(d), at this /, the corresponding tran- 
sition matrix elements are |toi| w 0.01, |<i2| ~ 0.07, 
and 1^02! ~ 0.13. When a microwave field is applied 
to drive a resonant transition |1) — > |2) = {£2), be- 



cause > » Fio at / ~ 0.493, the system can be 
pumped from |1) to |2) and then quickly decays to the 
ground state |0), while the process for thermally activat- 
ing the system from |0) to |1), via coupling to the environ- 
ment, will be very slow. Note that the coupling strength 
of the states |2) and |0) to the flux noise source is also 
proportional to the transition matrix element 1^20 1 j so the 
decay rate from |2) to |0) is proportional to (oc |i20p) 
according to the Fermi golden rule Therefore, the flux 
qubit is "cooled" because the population probability for 
the ground state |0) can be greatly increased, with re- 
spect to any unwanted excited state |1). Interestingly, 
this cooling mechanism corresponds to an "inverse pro- 
cess" of the usual state population inversion. For sim- 
plicity, here we use a weak microwave field. The driving 
field would need to be stronger to achieve cooling when 
the relevant transition matrix elements are small. This 
puts some constraints on the specific amplitudes used to 
achieve the desired result [l3|. Indeed, a recent experi- 
ment Q has successfully realized the microwave-induced 
cooling, lowering the temperature of a flux qubit relative 
to its surroundings. Thus, this microwave-induced cool- 
ing provides an efficient method for preparing the fiux 
qubit in its ground state. Below we use this prepared 
ground state to further cool a quantum system connected 
to the qubit. 

Cooling a quantum two-level system. — In the subspace 
spanned by |0) and |1), the fiux qubit (our AA) is mod- 
eled by Hq — ^huJioCTz- Here we consider a qubit-TLS 
system described by Ht = Hq + Htls + V + i?cnv, where 
Htls — \^^^'z is the Hamiltonian of a quantum TLS 
and i?env describes all the degrees of freedom in the en- 
vironment and their coupling to the TLS. Hereafter, the 
Pauli operators with primes refer to the neighboring TLS. 
The interaction Hamiltonian between the qubit and the 
TLS isV = %(cr+cr!_-hH.c.), with g = \tQi\/h. In the ex- 
perimental case (j), corresponding to Fig. [TJ^d), because 
1^01 1 is small at / 0.493, the coupling between the 
qubit and its environment is weak. To cool the TLS ef- 
fectively, after the qubit with fs = 0.22 is cooled to the 
ground state, we change the reduced magnetic flux fs to 
fs = 0.27, which corresponds to Fig^^lUe). For the qubit 
parameters used here, it is shown y that at / ~ 0.493 
(i.e., in between the level-crossing points), the adiabatic 
condition \h{Ei\{d/ dt)\Ej) / {Ei - Ej)\ < 1 can stifi be 
fulfilled for the three lowest levels by changing the ap- 
plied flux as fast as 0.1 $0 ns~^. This means that around 
this / the quantum states can be well preserved even 
when changing the flux very fast. More importantly, in 
the case of Fig.[TJe), because |toi| is much increased, then 
the qubit-TLS interaction Tig is strengthened by one or- 
der of magnitude. Here we assume that the quantum 
TLS is resonant to the qubit with fs = 0.27. Since the 
level spacing Tiwoi of the qubit with fs = 0.22 is different 
from that with fs = 0.27, thus at / ~ 0.493 the qubit 
with fs = 0.22 is off-resonant to the TLS. This gives an 
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FIG. 2: (Color online) (a) Ground-state probability pg ver- 
sus time t (in units of l/fa) and gr for Nt = Va/K, ~ 150. 
(b) Probability pg versus t and Nt for gr — 0.2tt. Here 
nth = 0.5 and Pe is chosen to be 0.4 at the initial time t — 0; 
g (r) is the interaction strength (time). 



even smaller effective qubit-TLS coupling. 

For simplicity, we now assume that the flux qubit is 
ideally cooled to the ground state |0) and then begins 
to resonantly interact with the quantum TLS at time U. 
When -ffonv is not included, the time evolution of the 
density operator of the TLS is governed by p(ti + t) = 
M{T)p{ti) and the gain operator is defined by M{t)p — 
Tr[exp(-iFT/?i)p(g) |0)(0| exp{iVT/n)], where Tr denotes 
the trace over the qubit states and r is the interaction 
time between the TLS and the flux qubit. 

When iJenv is considered, the dynamics of the density 
operator is described by 



dp 
'dt 



Ta ln[M(T)]p + Lp, 



(1) 



where is the rate for "switching-on" the AA to res- 
onantly interact with the TLS (each cycle includes the 
time required to cool the qubit) and L describes the dis- 
sipation of the TLS due to Hem- We model the environ- 
ment in Henv by a thermal bath. The operator L can be 
written as Lp = — ^^(nth + l)(o'+cr^p — cr'_pa'^) — 
\Knfti{a'_a'j^p— a'j^pa'J) -t- H.c, where k is the decay rate 
of the TLS and nth is the average number of bosons in the 
thermal bath (particularly, nth = at zero temperature) . 
Here we assume g > k, ensuring coherence between the 
qubit and its ancillary circuitry. 

For the neighboring quantum TLS, Eq. ^ can be ex- 
actly solved. The solution for pe = (e|p|e) is 



nth 
A 



exp(— A/it) 



nth 
A ' 



and Pg = {g\p\g) = 1 - Pe, where A = (2nth + 1) - 
TVt ln[cos^((7r)], with Nt — ra/n denoting the number of 
cycles for switching on the AA during the lifetime (= 
of the TLS. Because of the coupling to the qubit, 
the decay rate is now scaled by a factor A. Clearly, Pf, — 
nth/A and = 1 — nth/A at steady state. 

Because A is a periodic function of gT, both pe and 
Pg are also periodic; e.g., at gr = (2n — l)7r/2, with 
n = 1,2,..., A — > +CXD and Pg abruptly changes to 
Pg — 1; at gr — mr, with n = 0, 1, . . ., A = 2nth + l and pg 
slowly approaches = 1 — nth/(2nth + l). These features 
are clearly shown in Fig. [U^a) for pg with A^t = 150. To 
implement an efficient cooling, a smaller r is desirable, so 
we can only focus on the region gr £ [0, 7r/2] . Figure[2^b) 
shows the time evolution of Pg as a function of Nt for 
gr — 0.27r. Though gr is away from gr = 7r/2, one can 
still drastically cool the TLS by evolving Pg{t) to Pg ^ 1 
with a large Nt- 

Cooling a quantum resonator. — When the system con- 
nected to the flux qubit is a quantum resonator, the to- 
tal Hamiltonian becomes Ht = Hq + Hrcs + V + i/cnv, 
where iJrcs = huia^a describes the quantum resonator, 
V — ~hg{a+a -\- H.c.) is the interaction between them, 
and Hcnv describes all the degrees of freedom in the en- 
vironment and their coupling to the quantum resonator. 
Also, we assume that when cooling the quantum res- 
onator, the flux qubit is tuned in resonance to it. 

For the quantum resonator coupled to the flux qubit 
as well as to a thermal bath, the dynamics of the 
density operator of the quantum resonator is also de- 
scribed by Eq. ([1]). The operator L describes the dissi- 
pation of the quantum resonator induced by the thermal 
bath [l3 l: Lp ~ — ^^(nth + l){a'^ap + pa'^ a — 2apat) — 



^Knth(aaV + paa^ ^ 2a^ pa), where k is the damping 
rate of the quantum resonator and nth is the average 
number of bosons in the thermal bath coupled to the 
quantum resonator. In the present case, Eq. ([1]) can 
only be solved approximately. Here we use ln[M(T)] « 
[M — 1) — \{M — 1)^, which corresponds to neglecting 
terms of order 0{suY^{gTy/n)). The equation of motion 
for the boson number distribution p„ = (n|p|n) of the 
quantum resonator becomes 



dpn 
dt 



O-n+lPn+l — bn+lPn+2 — Cn+lPn 
-a-nPn + bnPn+1 + C„p„_l, 



(3) 



with a„ = raS{n)[l + ^S{n)] + K(nth + l)n, 6„ = 
^S{n)S{n + 1), and c„ — nrithn, where S{n) = 
sin^ {gr ^/n) . 

At steady state, dpn/dt — 0, which leads to a recursion 
relation for the steady boson number distribution p„: 



Pn-l = Pn 



(2) 



nth + 1 ^ NtS{n)[2 + S{n)] 
nth 2nnth 
At5(n)5(n + 1) 



-Pn+i- 



2nnth 



(4) 
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FIG. 3: (Color online) (a) Vacuum-state probability po versus 
time t (in units of 1/ra) and A'^t = Va/n. (b) Average boson 
number (n) versus t and Nt. Here nth = 0.5, gr = 0.27r, and 
(n) is chosen to be 1.5 at the initial time t = 0. 



where Nt — Va/ k, represents the number of cycles 
for switching on the AA during the lifetime of the 
quantum resonator. For TV ^ 1, pAr(nth + 1) S> 
PN+iNt S{N)S{N + 1)/2N, so we approximately have 

m-i = {^+ This is the initial 

condition for Eq. ^ and pat is determined by '^n^^Pn — 
1. 

Figure [3] displays the time evolution of both vacuum- 
state probability pq and average boson number (n) as a 
function of Nt for the quantum resonator, where gr = 
0.27r. As shown in Fig. ^a), one can evolve po{t) to 
Po ~ 1 using a large Nt- Figure [3Jb) shows that {n) ^ 
when Po ~ 1) revealing that the quantum resonator can 
also be effectively cooled. More interestingly. Figs. [D^b) 
and[2{a) give quite similar results, although very different 
models are used for the quantum systems connected to 
the qubit. This reveals that the cooling is robust and 
effective, irrespective of the chosen neighboring quantum 
systems. 

Discussion and conclusion. — The cooling approach 
studied here has potentially wide applications. For in- 
stance, the environmental noise is sometimes explained 
as mainly due to two-level fluctuators, in which one or 
a few fluctuators play a dominant role. Also, the en- 
vironment is often modeled by a boson bath, in which 
the bosons in resonance to the qubit play a dominant 
role. Here the quantum TLS can be used to model a 
two-level fluctuator and the quantum resonator can be 
used to model the dominant bosons of the environment 
in resonance to the qubit. Actually, the TLS defect that 
is most strongly coupled to the qubit may be off reso- 



nant to the qubit. If the off-resonance is large, the effect 
of the TLS on the qubit is not important. Otherwise, 
to cool the TLS defect, one can vary the reduced flux 
/ to tune the qubit to be in resonance with the defect. 
Also, when the environmental bosons in resonance with 
the qubit are cooled, one can tune the qubit by changing 
/ to further cool the off-resonant bosons. After cooling 
the dominant noise sources of the qubit, the quantum 
coherence of the cooled flux qubit will be enhanced. The 
quantum TLS can also model a solid-state qubit and the 
approach can be used to describe cooling two coupled 
qubits. Naturally, the quantum resonator can model a 
mechanical resonator at the nanometer scale. The cool- 
ing of mechanical resonators is currently a popular topic 
and its study provides opportunities to observe the tran- 
sition bewteen classical and quantum behaviors of a me- 
chanical resonator 13[. In our proposal, the quantum 



states can be manipulated quickly, due to the advantages 
of the proposed solid-state three-level system. Moreover, 
the cooling of both the flux qubit and the mechanical 
resonator can simultaneously enhance the quantum be- 
haviors of the two subsystems. This will help observe 
the transition between classical and quantum behaviors 
of the mechanical resonator via measuring the quantum 
states of the qubit. 

In conclusion, we have proposed an approach to simul- 
taneously cool a flux qubit and its neighboring quantum 
system. In each cycle of cooling, the flux qubit is first 
prepared to the ground state, following an inverse process 
of the state population inversion, and then switched on 
to resonantly interact with the neighboring quantum sys- 
tem. As typical examples, we model the quantum system 
connected to the qubit by either a TLS or a resonator. 
Our results show that the cooling is robust and effective, 
irrespective of the chosen quantum systems. 
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